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QQ Abstract 
^. 

Direction measurement of weakly interacting massive particles in time-projection 
chambers can provide definite evidence of their existence and help to determine their 
properties. This article demonstrates several concepts for charge amplification in 
time-projection chambers that can be used in direction-sensitive dark matter search 
• 5^ experiments. We demonstrate reconstruction of the "head-tail" effect for nuclear 

^ recoils above 100 keV, and discuss the detector performance in the context of dark 

5^ matter detection and scaling to large detector volumes. 
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1 Introduction 



Weakly interacting massive particles (WIMPs) are good candidates for con- 
stituents of a dark-matter halo around our galaxy, and a target of several 
experimental searches (see e.g. [1] for review). An annual oscillation in the 
signal at the level of 1-2% was recently confirmed by DAMA collaboration [2] . 
If the asymmetry is due to Earth's motion around Sun then much larger diur- 
nal asymmetry in the WIMP direction is expected [3]. However, the DAMA 
and other experiments that report zero signal use liquid and solid detec- 
tor materials that allow for more compact detectors, but destroy information 
on the WIMP direction and the sense of direction ("head-tail"). Therefore, di- 
rectionality is particularly valuable, both for suppressing background and for 
confirming the DAMA result by correlating a candidate dark-matter signal 
with astrophysical phenomena. 

Reconstruction of the WIMP direction can be accomplished with detectors 
using low-pressure gas as the target material. Several groups, DRIFT [6], 
NEWAGE [7j, MIMAC |8] and DMTPC [9J have reconstructed low- momentum 
recoils created in elastic neutron scattering using low-pressure time-projection 
chambers for dark matter searches. The reconstruction of the sense of direc- 
tion ("head-tail") using the scintillation profile of a recoiling track has been 
only recently achieved [S]. 

The low density of gaseous detectors and the small WIMP cross section ne- 
cessitate the use of large detector volumes with fine detector granularities. 
A significant improvement of the current experimental limits may require a 
ton-scale detector. In the case of CF4 as the detector material, one ton of gas 
occupies a volume of approximately 16 x 16 x 16 m'^ at 50 Torr of pressure. A 
50 keV fluorine recoil created in a WIMP collision travels 1.5 mm at 50 Torr of 
pressure. Therefore, a multi-cubic meter detector with resolution of the order 
of hundreds of micrometers is needed for a directional dark matter experiment, 
requiring further progress in detector technology in order to observe WIMPs. 



2 Detector designs 

In this paper we propose and demonstrate several designs for time-projection 
chambers that can be used in large-volume, directional dark matter search 
experiments with "head-tail" discrimination. In a time projection chamber, 
a WIMP creates a nuclear recoil that makes electron-ion pairs as it slows 
down in the detector gas. We use CF4 gas that has good charge multipli- 
cation and scintillation properties piHlllll2j . Fluorine, in addition, has non- 
zero angular momentum, which allows probing for spin- dependent dark mat- 
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ter interactions. The electron diffusion limits the maximum drift distance to 
about 25 cm [TIP] . The initial ionization electrons drift in electric field to- 
ward the charge-amplification region, where in case of CF4 , the threshold 
for charge multiplication is approximately 45 V/(cm-Torr) [10]. Traditionally, 
design choices for amplification regions have been based on the multiwire pro- 
portional chamber (MWPC) [13], where strong electric fields are created in 
the vicinity of thin wires. The pitch between wires, and therefore the spa- 
tial resolution in MWPC is limited to above 1-2 mm due to mechanical and 
electrostatic reasons. Finer spatial resolution can be achieved with micropat- 
tern detectors (e.g. [BJ), but the size and the gain of detector modules are 
constrained by manufacturing limitations. 

In this paper we demonstrate designs for micropattern detectors that provide 
good gain, fine granularity and scalability to large volumes. Electrodes used in 
the charge multiplication are made of woven meshes, indium-tin-oxide (ITO) 
films, or copper. All designs allow production of large-area modules as the 
meshes and ITO foils are produced in 1.2 m-wide rolls. In the first case shown 
in Figure [1^), the amplification region is made of a stainless-steel mesh and 
copper-clad GIO board. The stainless steel mesh is made of 28 /im diameter 
wires with periodicity of 256 /im, which gives optical transmittance of 77%. 
Woven meshes are used in air filtration systems and can be mass produced 
in a cost-effective way. The mesh and the copper-clad board are separated 
by fluorocarbon resistive (fishing) wires of 0.54 mm diameter and spaced ev- 
ery 2 cm, with transparency of 97%. With this design, we demonstrate the 
reconstruction of 2D recoil segments with improved gain and no additional 
cost compared to our previous detector. Further improvement in gain may be 
achievable with the use of transparent electrodes that allow scintillation light 
to be read out from two sides. We test a design shown in Figure [T|d), where 
the amplification electrodes are made of ITO layers deposited on mylar foil 
and a stainless steel mesh. ITO is 90:10 mix ratio, by weight, of indium-oxide 
(In203) and tin-oxide (Sn02), respectively. ITO is widely used in consumer 
electronics (e.g. touch screens, LCD screens) and available in 1.2 m wide rolls of 
ITO-coated mylar foils. The thickness of the ITO coating determines electrical 
and optical properties of the foil: we use a surface resistance of 15 Ohm/cm^, 
which gives light transmittance of around 80%. In the third design shown in 
Figure [T]^), we use two meshes to create the amplification field and the same 
fluorocarbon wires to separate the planes. The transparency of electrodes can 
potentially increase the gain by almost a factor of two due to simultaneous 
readout of two drift volumes. 

The experimental setup with a time-projection chamber (TPC) and a CCD 
camera for optical readout [15] of the amplification plane is shown in Figure |2} 
The cathode mesh used to create the drift field has a periodicity of 312 fim 
and wire thickness of 31 /xm, which gives an optical transparency of 88%. The 
drift distance is limited to approximately ±5 cm by the size of the vacuum 
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vessel. The CCD camera is manufactured by Finger Lake Instrumentation and 
equipped with a Kodak KAF-0401ME chip with a cooler that maintains the 
temperature in the range [—20,-18] C to minimize electronic noise (25 ADU). 
The photographic lens has an aperture ratio, f/# of 0.95, and a focal length 
of 25 mm. 
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Fig. 1. Schematics of amplification regions formed with a) mesh and a copper sheet, 
b) mesh and ITO-film, and c) two meshes. In all cases the electrodes are separated 
with 0.54 mm fluorocarbon wires. 




Fig. 2. A schematic of the TPC chamber used in these studies. The inner diameter 
of the chamber is 23 cm, and the height of a drift region is 5 cm. The ITO-mesh 
configuration has two drift regions that can be read out simultaneously. 
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3 Results 



3.1 Detector gain 

We measure the detector gain using 5.5 MeV a particles from a collimated 
^^^Am source. The chamber is filled with CF4 at various pressures in the range 
50-200 Torr. The drift field is 500 V/cm, and the amplification voltage ranges 
between 0.65 and 1.1 kV. Images are taken sequentially with 500 ms expo- 
sure time and resolution of 96x64 pixels, where each pixel is 488 x 488 /im^ 
(72x72 /im^ on chip). Pixels that have intensity greater than 5 standard de- 
viations from the mean dark field at least 10% of the time are fiagged as 'hot 
channels' and excluded in the data analysis. ADC bias is corrected for by 
subtracting from each image the average of 100 images taken with the shutter 
closed. 

The gain of the detector is determined from the intensity of scintillation light 
recorded by the CCD camera in a 5 mm track segment close to the alpha 
source. The stopping power of alpha tracks in this region is approximately uni- 
form and we estimate it using the SRIM [17] program, as 0.75 MeV/ (mg/cm^), 
which corresponds to 87, 130, 175, 264 and 355 keV of energy loss at 50, 75, 
100, 150 and 200 Torr of CF4 pressure, respectively. 

We convert the observed light intensity into charge gain by accounting for 
the transmittance of the amplification and drift meshes (70%), vessel window 
(90%), camera lens and window (90%). Scintillation light created in the lower 
drift region has an additional transmittance loss of approximately 60% as it 
passes through ITO layers and acrylic plate, or two amplification meshes. We 
compute the lens acceptance to be 1.1 ■ 10"'^ and estimate the average CCD 
efficiency as 40 ± 10% using the manufacturer's quantum efficiency curve and 
the scintillation spectrum of CF4 gas [12]. The gain of the camera is measured 
to be 1.6 ADU/e~. The ratio of the total number of scintillation photons to 
the electrons in the avalanche is taken as 1/3 [TTP^ . We use w = 54 eV for 
the average ionization energy in CF4 gas, or 18.5 ionization e~ per keV of 
energy loss [16]. Hence, the total gain, g is computed from observed number 
of counts per keV of energy loss, Iccd as g = Iccd/'^-^ ■ 10~^. 

We measure the gain in the mesh-copper detector (Figure [T^) and plot the 
observed intensity Iccd as a function of the amplification voltage and pressure 
in Figure [3^). The voltage is increased until the total charge created in the 
amplification region reaches a sparking threshold. The maximum gain varies 
from 15-45 ADU/keV, which is roughly 13 times larger than charge amplifi- 
cation achieved with wires [5], after accounting for different light collection 
efficiencies. Since the lower ionization density allows larger charge multiplica- 
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tion, the maximum gain increases with decreasing pressure, as is evident in 
Figure [3^). The maximum electronic stopping power of alpha particles occurs 
in the Bragg peak at around 800 keV. The same electronic stopping power is 
reached at 250 keV for fluorine ions, presenting an upper-energy threshold for 
detection of nuclear recoils at a given gain. We also compare different gaps 
between anode and ground planes using different widths of fishing lines. After 
adjusting the amplification voltage, we find that the separation has no effect 
on the maximum achievable gain, which is limited by the total charge created 
in the amplification region. 



The gain (or reduction in the number of CCD cameras) can be further im- 
proved by placing drift regions opposite to each other along with transparent 
amplification electrodes, as shown in designs with ITO and all-mesh electrodes 
in Figures Qi,c). Two Am sources are placed in opposite drift regions ap- 
proximately 1 cm from the amplification region, which is made of an ITO-mesh 
sandwich (Figure [T]d). Both drift regions are read out by the same CCD cam- 
era as shown in Figure [2j An image taken with a 500 ms-exposure of several 
alpha tracks from two ^^^Am sources is shown in Figure m Note that the sig- 
nal from tracks in the lower source is attenuated due to passage through two 
ITO-coated anodes and a ground mesh. We measure the attenuation to be 
(52 ± 2)%, which is close to the expected attenuation of (58 ± 1)% when we 
assume the same gain for the upper and lower amplification regions. We eval- 
uate the gain using the signal from the upper source, and plot it in Figure [sjo). 
The gain is found to be 20-30% lower than with the copper-mesh amplification 
plane. Further improvements are possible with foils that have higher resistivity 
(light transmittance) and better quality control during shipment and module 
production, as our current foils arrived slightly damaged. 



Finally, we test a detector that has the amplification plane made entirely of 
meshes. The gain measurement is shown in Figure [sj:;) . Values are comparable 
with ITO-based detector and approximately 30% lower than the copper-mesh 
detector. Operation with a readout of two drift regions can be realized by 
adding a third mesh. 



An energy resolution of approximately 10% and spatial resolution of 400 fim 
are measured here and found to be comparable with our previous measure- 
ments using MWPC |9]. Contribution to the spatial resolution from the dif- 
fusion of electrons that are 2 cm above the amplification plane is 140 /im, 
finite CCD bin size adds 488/-\/T2 ^ 140 /xm, and finite mesh pitch adds 
256/v^ ^ 74/im. 
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Fig. 3. The gain measured in the three amphhcation designs. 
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3.2 Gain uniformity 



The uniformity of the gain is demonstrated by observing a clear Bragg peak 
in mesh-ITO and mesh-copper detector designs (Figure [l^,b). A scintillation 
profile of an alpha track that has a range of approximately 4.5 cm in 250 Torr 
of CF4 is shown in Figure |4] 

In the case of the mesh-mesh detector (Figure [Tj^), gain non-uniformity can 
occur on a smaller scale if the relative alignment between the meshes varies 
across the surface of the detector. The charge multiplication and light collec- 
tion are more efficient if wires of the anode mesh fall in between wires of the 
ground mesh, i.e. with offset of half of a mesh pitch in both x and y direc- 
tions. The variation in relative alignment between the two meshes may be the 
result of a non-zero relative angle, non-uniform mesh tensioning, or variation 
in spacing between the mesh planes. We show an example of this by taking a 
photograph of light scattered from the bottom mesh and transmitted through 
the upper mesh, as shown in Figure |5^). A moire pattern with periodicity 
of roughly 2 mm is created by the change in relative alignment between two 
meshes. In order to evaluate the effect on the gain, we expose the detector to 
x-rays from an ^^Fe with gas pressure set to 400 Torr and anode voltage to 
1.74 kV. The average gain of approximately 9 ADU/keV is determined from 
the total light emitted by 5.5 MeV alpha particles. An ^^Fe source is mounted 
on the top of the cathode mesh and several minutes of 1-second exposures 
are taken, with accumulated CCD image shown in Figure |5|3). The average 
light intensity decreases with distance from the source due to the reduced 
flux of x-rays. The semi-circular shadow on the right is due to the Fe source. 
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Fig. 4. An image of alpha tracks from 500 ms exposure at 200 Ton with two Am 
sources placed in opposite drift regions. The amphfication electrodes are mesh and 
ITO foil. 
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and the long shadows are due to fishing fine separators. The checkered hght 
pattern seen in the plot is due to variation of the gain and follows the moire 
pattern. Using pixels that have equal distance from the source, we estimate 
a gain variation to be approximately 50% between maxima and minima. The 
non-uniformity in gain can be minimized, if not avoided, by applying extra 
care during fabrication of detector. Elimination of the moire pattern during 
the mesh stretching and assembly can be used as a monitor for quality control. 
Any remaining gain variation can be accounted for by calibrating with x-ray 
sources. 



a) Moire pattern b) ^^Fe exposure 




Fig. 5. Evaluation of gain uniformity in the mesh-mesh detector using moire pattern 
and ^^Fe exposure. The distance between horizontal resistive wires is 2 cm. 

3.3 "Head-tail" effect 

Finally, we demonstrate the ability to determine the directional sense of low- 
momentum nuclei, the "head-tail" effect. The ionization rate by a low-momentum 
nucleus decreases as it slows down in the detector gas, and the sense of its di- 
rection can be deduced from the scintillation profile along the track. We place a 
^^^Cf source approximately 2 m from the detector's view-field of 4.8 x 3.2 cm^, 
resulting in a well determined direction of the neutron fiux. We use the copper- 
mesh amplification plane (see Figure [l^) with the anode voltage set to 740 V 
and the CF4 pressure at 75 Torr. Typical ranges for nuclear recoils created in 
neutron scattering are of the order of a few millimeters. 

Taking 6000 1-second exposures without a trigger, we search for clusters of pix- 
els that are 3 standard deviations above the CCD noise. In most cases there is 
only one cluster per event, but otherwise we select the most energetic cluster. 
The energy is computed from the sum of pixel yields after background subtrac- 
tion, using calibration resulting from inserting an alpha source three times: 
before, in the middle and after the neutron exposure. We find an average gain 
of 27 ADU/keV. The gain decreases by 8% during the measurements, probably 
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due to air leakage into the vacuum vessel. We describe quenching correction 
for a recoil with energy E in terms of electronic Se{E) and nuclear Sn{E) 
stopping powers, q{E) ^ (^e(^) + ^n(^) ■ 0.3) / (^e(^) + 5'„(E)), where fac- 
tor 0.3 is chosen to be small and non-zero. The quenching correction, q, which 
is about 25% at 100 keV and 3% at 900 keV [T7|IT8] , is used in the calculation 
of light output for fluorine recoils. The range of the recoil track is determined 
as the maximum distance between the pixels in the cluster. Recoils are de- 
fined to have at least 6 CCD bins that are not touching the boundary of the 
CCD view field. We require a ratio of principal moments of inertia greater 
than 5, which imposes a minimum energy cut of around 100 keV on nuclear 
recoils. Discharge events are removed with a cut on the maximum light per 
event, and remaining spot-like events are removed with the cut on the ra- 
tion of principal moments. A plot of energy versus range of simulated events 
and data is shown in Figure |6] for recoil candidates that have reconstructed 
energy above 100 keV. The simulation is based on cross sections found in a 
nuclear scattering library [19] and fluorine ion propagation determined from 
the SRIM [IT]. 
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Fig. 6. Range (top left), skewness (top right), and cosine of 2D recoil angle (bottom 
left) vs. reconstructed energy for nuclear recoil candidates in a ^^^Cf exposure at 
75 Torr. Black points are data, the box-histogram is simulation. Signed distribution 
of the cosine of 2D recoil angle (bottom right), with data and simulation normalized 
to the same area. 

We show images of recoil tracks with highest energies in Figure [7j In all images 
neutrons are incident from the right, along the x axis. The nuclear recoils also 
propagate from the right, and we observe a decreasing light intensity as the 
recoils slow down in the CF4 gas, as expected. We quantify the scintillation 



10 



asymmetry along the track using skewness, which is described previously in 
more detail [U]. It is a dimensionless ratio of the third moment and the root- 
mean-square cubed of the light distribution along the track. Negative skewness 
is expected for recoils heading in the direction of the neutron beam. The top- 
right plot in Figure |6] shows the skewness distribution for tracks in data and 
simulation. 

The two pixels with maximum separation that are used in the range measure- 
ment are also used in the calculation of the recoil angle. Figure [8] shows the 
resolution for the recoil angle based on the simulation. The sense of direction 
is determined from the "head-tail" (skewness) measurement. We show for the 
first time, distributions for the signed cosine of the 2D recoil angle, cos OnecoU 
vs. energy and signed cosine of the 2D recoil in Figure |6| 

We define a quality factor, Qht for the "head-tail" asymmetry at given recoil 
energy, Efi as 

Qht{Er) = e{ER) ■ (1 - 2 uj{ER)f (1) 

where the recoil reconstruction efficiency, e and the fraction of wrong "head- 
tail" assignments, uj are determined from simulation. The Qht is the effective 
fraction of reconstructed recoils with "head-tail" information, and the error 
on the "head-tail" asymmetry scales as 1/ \fQuT- Figure [s] shows the quality 
factor in the energy range that has been explored so far. 



4 Conclusion and outlook 

Directional detection of dark matter requires large detector volumes with fine 
granularities. TPC modules filled with low-pressure CF4 gas can provide di- 
rectionality and "head-tail" directional sense, in particular for spin-dependent 
dark matter searches. We have demonstrated three possibilities for charge- 
amplification that allow directional detection and "head-tail" determination. 
In all three cases the gain is improved by more than an order of magnitude 
compared to the previous MWPC design [9]. This results in improvement 
in the "head-tail" discrimination due to the larger gain, 2D detection of re- 
coil images, and reduced pressure from 200 to 75 Torr, allowing for longer 
recoil tracks. We estimate the sensitivity for the WIMP detection using stan- 
dard assumptions about the dark matter halo [2T]. Cross-section limits for 
spin-dependent WIMP scattering on proton are shown in Figure |9j We can 
improve current experimental limits |1] with approximately 0.1 kg • y of CF4 
exposure, and test MSSM models [22] with approximately 100 kg ■ y. In both 
cases we assumed operation in an underground laboratory with a neutron- 
induced background rate of 0.01 events/(keV • kg • y) [20], uncorrelated with 
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Cygnus direction, and a 50 keV recoil-energy threshold. We plan to pursue 
these technologies and construct a cubic-meter module that will be a basic 
building block of a ton-scale detector. 
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